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Numerical study of order parameter dynamics in 
the course of second order (Landau-Ginzburg) symmetry 
breaking transitions shows that the density of topological 
defects, kinks, is proportional to the fourth root of the rate 
of the quench. This confirms the more general theory of 
domain-size evolution in the course of symmetry breaking 
transformations proposed by one of us Using these 

ideas, it is possible to compute the density of topological 
defects from the quench timescale and from the equilibrium 
scaling of the correlation length and relaxation time near 
the critical point. 
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The dynamics of symmetry breaking phase transi- 
tions has been the focus of attention because of its im- 
plications for cosmological scenarios and its 
importance in the context of condensed matter physics 
|,§|J|,|l^,|ri|. The likely cosmological setting for this 
process involves a second order phase transformation 
associated with a non-conserved order parameter (p. 
Initially, the system is near equilibrium in the symmet- 
ric phase with {ip) = 0. When the effective potential 
assumes the "sombrero" shape, (p is forced to make a 
choice of one of many possible broken symmetry true 
vacua. As pointed out by Kibble |^ in the cosmological 
context, the finiteness of the speed of light implies that 
the choices of the new broken symmetry vacuum must 
be independent in sufficiently distant regions of space. 
Therefore, when the homotopy group, characterizing 
the relation between the manifold of broken symmetry 
vacua and the space in which tp evolves, is non-trivial, 
locally independent choices of the new vacuum will lead 
to topologically stable objects such as monopoles, cos- 
mic strings, or domain walls. The initial density of 
such defects is of great interest and cannot be deduced 
from the scenario sketched above, save for a much too 
generous lower limit which follows from the lightcone 
causality alone | |T^ ]. 

As suggested in [Q , it is also possible to study "cos- 
mological" mechanisms of defect formation in con- 
densed matter. In this context, the estimate of defect 
density based on speed of light arguments is consider- 
ably less useful than in cosmological situations. Thus, 
along with the proposal for condensed matter tests of 
defect formation in cosmological scenarios, a theory 



aimed at computing the defect density was developed 
[Pjlljl . The key idea is to realize that the order param- 
eter can only become correlated through its dynamics. 
Furthermore, in the vicinity of the critical temperature 
Tc, order parameters exhibit universal behavior char- 
acterized by two simultaneously occurring divergences. 
That is, when parametrized in terms of the relative 
temperature e — (Tc — T)/Tc, both the equilibrium 
healing length (also known as the correlation length) 



and the dynamical relaxation time 
r = ro/\er 



(1) 



(2) 



will simultaneously diverge at the critical temperature 
(see Fig. p. Above, and tq characterize the low 
temperature (T = 0, e = 1) values of the healing length 
and relaxation time, respectively. 

The healing length and the dynamical relaxation 
time have similar physical significance. ^ is the dis- 
tance over which the order parameter returns to its 
equilibrium value when perturbed, for instance, by the 
boundary conditions. ^ is also the typical scale of the 
perturbations; that is, to the leading order, the corre- 
lation function of perturbations away from equilibrium 
value, 6ip{x,t), behaves as 



{Sipix,t),Scp{x + A,t)) ^ exp(-|A|/e) 



(3) 



Below Tc, the spatial order is established on scales 
much larger than ^; however, the scale over which the 
healing occurs — for example, from the "wounds" in- 
flicted by topological defects — is characterized by ^. 
Similarly, r characterizes the time required for the or- 
der parameter to relax to its equilibrium value. Dur- 
ing the phase transition, in the immediate vicinity of 
the critical temperature, the motion will be often over- 
damped (that is, dominated by the first time derivative 
of the order parameter) and r = To/|e|. 

The estimate of the density of defects put forward in 
Ref. ||] is based on a linear quench. 



t/TQ 



(4) 



which is expected to be a suitable approximation in the 
neighborhood of Tc- In (|^), tq is the quench timescale 
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and t the time before {t < 0) and after {t > 0) the 
transition (e = 0). When the critical temperature is 
crossed at the fixed rate given by Eq. there will 
come a moment when the order parameter — because 
of the critical slowing down implied by Eq. — will 
be simply unable to adjust its actual correlation length 
to the equilibrium ^ given by Eq. (|). This wiU occur 
when the time remaining until the transition, i, equals 
the relevant dynamical relaxation time T{t), that is, 
when 

r(t) = To/{\i\/TQr = i . (5) 

This equation yields a relative temperature 

e = i/TQ = {To/TQ)^, (6) 

at which the evolution of the order parameter will cease 
to follow equilibrium. At this point, the correlation 
length of (fi will cease to diverge in accord with Eq. (|^) , 
as the phase transition region is traversed. Instead, ^ 
will reach a value approximately given by 

e = eo/|er =Co(rQ/To)T*7r . (7) 

Regions more distant than ^ will be forced to select the 
new vacuum independently. 

An essentially equivalent estimate of the character- 
istic domain size is obtained from the "causal horizon" 
idea, but only when the appropriate velocity of the sig- 
nal propagation on ip is adopted [0,0 . That velocity is 
of the order of w = ^/r for the dynamical processes cor- 
relating different regions of the order parameter. Thus, 
the size of the corresponding causal horizon size (and, 
therefore, the size of the domain) will be no more than 
h{i) = v{i)i ~ |. 

In any case, the size of independently selected do- 
mains of the new vacuum will be given, to the leading 
order, by ^, which will also determine the typical sep- 
aration of topological defects and, therefore, the ini- 
tial defect density n. For instance, for monopoles n « 
l/{fO^ , with D the space dimension and / ^ 0{1) 
a factor presumably somewhat larger than unity. / 
takes into account the possibility that the independent 
choices of the vacuum may be similar anyway, and that 
the correlations length will slowly grow following the 
instant t due to diffusion, etc.; the key prediction of 
the theory of 0] is in any case the scaling, as is usually 
the case in the critical phenomena. 

The aim of the work in this Letter is to test this 
theory with a computer experiment. The immediate 
motivation of our research comes from recent super- 
fluid experiments [f7|,p|jTc|] . These experiments follow 
the original suggestion in Ref. [Q and appear to sup- 
port the estimate of defect formation based on ^. While 
the quench-generated density of vortex lines is some- 
what uncertain, it is, nonetheless, in accord with the 
theory summed up above but in conflict with the ideas 



based on activation and Ginzburg temperature p^ . 
However, these experiments were, at least so far, un- 
able to vary the quench time tq. Thus, it was not 
possible to systematically test the key predictions of 
Ref. namely the power law dependence of the size 
of the fragments of the broken symmetry vacuum on 
the quench rate Tq^ given by Eq. (0), and the comple- 
mentary dependence of the initial number of defects. 

To investigate this issue, we considered the numerical 
evolution of a one-dimensional system for a real 
field ip according to the equation of motion derived 
from the Landau-Ginzburg potential V{ip) = {ip'^ — 
2eip'^ + l)/8 . The system is assumed to be in contact 
with a thermal reservoir. Thus, it obeys the Langevin 
equation 

<p + r](p- dj;.j;ip + d^V{ip) = 1? . (8) 

The noise term has correlation properties 

(t9(a;, t),^{x',t')) ^ 2ri0dix' - x)5{t' - t) . (9) 

In Eq. (^, 77 is the overall damping constant which also 
helps characterize the amplitude of the noise through 
Eq. (^). The parameter e measures the distance from 
the phase transition and varies according to Eq. (^; 
that is, e = min(l,i/rQ). 9 describes the temperature 
of the reservoir and is kept constant throughout this 
work. This separate parametrization of e and 9 was 
adopted to correspond to the situation in superfluid 
He^, where the symmetry breaking can be induced by 
the change of the pressure and occurs with inconse- 
quential adjustments of the absolute temperature 0,0]. 
The second order time evolution allows us in principle 
to make contact with cosmology, but, in the regime 
considered here, the evolution is dominated by the dis- 
sipative term rjip. Hence, in effect, we are dealing with 
the time dependent Landau-Ginzburg equation. 

We investigate the creation of kinks as a function 
of Tq by starting at some e < suitably above the 
transition, and then gradually adjusting the shape of 
the potential in accord with Eq. (Q). Figure || shows 
a sequence of "snapshots" of (p obtained in the course 
of such a quench. When e < 0, (/? fluctuates around 
its expectation value = 0. The same situation ini- 
tially persists for slightly positive e. However, further 
below Tc , ip settles locally around one of the two alter- 
natives: {ip) = ■ Moreover, local choices of one of 
these two alternatives cannot be easily undone once a 
certain symmetry breaking is selected, unless the kinks 
are separated by distances no larger than ^ (see Fig. 

To test the predictions of defect density in Ref. M, 
we note again that, for a sizable r/ and sufficiently 
small e, the damping term 7](p in Eq. (H) is bound to 
dominate. Such overdamped evolution will take place 
whenever rj^ > \e\ [Q. In our case, 77 > 1 and 
|e| ^ 1. The characteristic relaxation time is then 
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given by r — ri/\e\ . Consequently, fi — I in Eq. (|^), 
and 

i=Vv^- (10) 

The corresponding value of the relative temperature 
immediately yields, e = t/xQ = y^ij/rQ . It follows 
that 

e = Co(rQ^)l/^ (11) 

where we have adopted — 1/2 in accord with the 
Landau-Ginzburg theory and in agreement with the 
critical exponents inferred from the behavior of the 
healing length in Fig. |. Eqs. @ and (|ll|) are ex- 
pected to be applicable as long as the condition rj^ > |e| 
holds at e, which in turn implies jy/rg > 1. 

Figure ^ is the principal result of our paper. It il- 
lustrates the density of kinks obtained in a sequence 
of quenches with 77 = 1 for various tq values. For 
each Tq, the phase transition was simulated 15 times, 
starting at e = — 1 (except for the shortest and longest 
TQ, which were initiated at e = — lO/y'fg). Our com- 
putational domain had periodic boundary conditions. 
Production runs were carried out with resolution of 
16,384 gridpoints. Convergence was checked by com- 
paring results at different resolutions. The "physical" 
size of the computational ring in the production runs 
was 2,048 units. At this scale, the ring is large enough, 
so boundary effects are avoided. 

The number of kinks produced by the quench was 
obtained by counting the number of zeros of if. Above 
and immediately below Tc , there is a significant num- 
ber of zeros which have little to do with the kinks (see 
Fig. 1^). However, as the quench proceeds, the number 
of zeros quickly evolves towards an "asymptotic" value 
(see Fig. ||). This change of the density of zeros and 
its eventual stabilization is associated with the obvious 
change of the character of (p and with the appearance 
of the clearly defined kinks. By then, the number of 
kinks is nearly constant in the runs with long tq, and, 
correspondingly, their kink density is low. Even in the 
runs with the smallest tq, there is still a clear break 
between the post-quench rates of the disappearance of 
zeros and the long-time, relatively small rate at which 
the kinks annihilate. 

In the regime here investigated, the theoretical scal- 
ing relation for the number density of kinks, n « 
(v/'^q)^ ~ I appears to be well satisfied, with 

/ ~ 8. It is interesting to note that a similar value 
of / has been deduced from the experiment in Ref. 
0. Experimentally we find n = (8.7 x 10"^ ± 1.0 x 
10"^)Tg°-^*^""^, when the kinks are counted at ap- 
proximately the same I/tq value. A similar scaling is 
also obtained for kinks counted at equal t times after 
the quench. 

In summary, our numerical experiment appears to 
provide a strong confirmation of the theoretical pre- 



dictions given by one of us [|l|. As expected, the den- 
sity of topological defects is somewhat less than the 
inverse of ^ but of the right order of magnitude. Most 
importantly, the scaling of the kink density with tq 
follows closely theoretical expectations. These results 
are also supported by the dependence of the number 
of kinks on the damping parameter r/, as well as by 
the preliminary results of the computer experiments 
involving complex order parameter and/or more than 
one dimensional space |l^ . 
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FIG. 3. Number of defects as a function of quench 
timescale. The straight line is the best fit to n oc Tq" 
with a — 0.28 ± 0.02 (x^ ~ 1.96). This exponent compares 
favorably with the theoretical prediction of 1/4 based on 
the theory in Ref. Q. 



FIG. 1. Characteristic equilibrium correlation (or heal- 
ing) length, ^ = Co/|er, for the system under investigation. 
The best fitting yields = 1.38 ± 0.06, f = 0.41 ± 0.03 
(X^ = 1.2) above Tc, and Co = 1-02 ± 0.04, u = 0.48 ± 0.02 
(X^ = 3.7) below Tc, close to the Landau- Ginzburg expo- 
nent of v = 1 /2. 
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FIG. 2. Snapshots of ip during kink formation with a 
quench timescale of tq — 64 and damping parameter 
rj = 1. The figures, from top to bottom, correspond to 
t = -80.0, 7.5, 32.5, 45.0 and 333.0, respectively. 




FIG. 4. Average number of zeros as a function of time in 
units of tq; from top to bottom tq = 4, 8, 2048, 4096. 
The number of kinks used in Fig. ^ were obtained at 
t/TQ ~ 4, except in the large, computationally expensive, 
TQ cases where an extrapolated value was used. 
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